Disregarding multimappers leads to biases in the functional assessment of NGS data.

BACKGROUND: Standard ChIP-seq and RNA-seq processing pipelines typically disregard sequencing reads whose origin is ambiguous ("multimappers"). This usual practice has potentially important consequences for the functional interpretation of the data: genomic elements belonging to clusters composed of highly similar members are left unexplored. RESULTS: In particular, disregarding multimappers leads to the underrepresentation in epigenetic studies of recently active transposable elements, such as AluYa5, L1HS and SVAs. Furthermore, this common strategy also has implications for transcriptomic analysis: members of repetitive gene families, such the ones including major histocompatibility complex (MHC) class I and II genes, are under-quantified. CONCLUSION: Revealing inherent biases that permeate routine tasks such as functional enrichment analysis, our results underscore the urgency of broadly adopting multimapper-aware bioinformatic pipelines -currently restricted to specific contexts or communities- to ensure the reliability of genomic and transcriptomic studies.

Three concurrent mechanisms generate gene copy number variation and transient antibiotic heteroresistance.

Heteroresistance is a medically relevant phenotype where small antibiotic-resistant subpopulations coexist within predominantly susceptible bacterial populations. Heteroresistance reduces treatment efficacy across diverse bacterial species and antibiotic classes, yet its genetic and physiological mechanisms remain poorly understood. Here, we investigated a multi-resistant Klebsiella pneumoniae isolate and identified three primary drivers of gene dosage-dependent heteroresistance for several antibiotic classes: tandem amplification, increased plasmid copy number, and transposition of resistance genes onto cryptic plasmids. All three mechanisms imposed fitness costs and were genetically unstable, leading to fast reversion to susceptibility in the absence of antibiotics. We used a mouse gut colonization model to show that heteroresistance due to elevated resistance-gene dosage can result in antibiotic treatment failures. Importantly, we observed that the three mechanisms are prevalent among Escherichia coli bloodstream isolates. Our findings underscore the necessity for treatment strategies that address the complex interplay between plasmids, resistance cassettes, and transposons in bacterial populations.

Integrative and Conjugative Elements and Prophage DNA as Carriers of Resistance Genes in Erysipelothrix rhusiopathiae Strains from Domestic Geese in Poland.

Goose erysipelas is a serious problem in waterfowl breeding in Poland. However, knowledge of the characteristics of Erysipelothrix rhusiopathiae strains causing this disease is limited. In this study, the antimicrobial susceptibility and serotypes of four E. rhusiopathiae strains from domestic geese were determined, and their whole-genome sequences (WGSs) were analyzed to detect resistance genes, integrative and conjugative elements (ICEs), and prophage DNA. Sequence type and the presence of resistance genes and transposons were compared with 363 publicly available E. rhusiopathiae strains, as well as 13 strains of other Erysipelothrix species. Four strains tested represented serotypes 2 and 5 and the MLST groups ST 4, 32, 242, and 243. Their assembled circular genomes ranged from 1.8 to 1.9 kb with a GC content of 36-37%; a small plasmid was detected in strain 1023. Strains 1023 and 267 were multidrug-resistant. The resistance genes detected in the genome of strain 1023 were erm47, tetM, and lsaE-lnuB-ant(6)-Ia-spw cluster, while strain 267 contained the tetM and ermB genes. Mutations in the gyrA gene were detected in both strains. The tetM gene was embedded in a Tn916-like transposon, which in strain 1023, together with the other resistance genes, was located on a large integrative and conjugative-like element of 130 kb designated as ICEEr1023. A minor integrative element of 74 kb was identified in strain 1012 (ICEEr1012). This work contributes to knowledge about the characteristics of E. rhusiopathiae bacteria and, for the first time, reveals the occurrence of erm47 and ermB resistance genes in strains of this species. Phage infection appears to be responsible for the introduction of the ermB gene into the genome of strain 267, while ICEs most likely play a key role in the spread of the other resistance genes identified in E. rhusiopathiae.

Spatially revealed roles for lncRNAs in Drosophila spermatogenesis, Y chromosome function and evolution.

Unlike coding genes, the number of lncRNA genes in organism genomes is relatively proportional to organism complexity. From plants to humans, the tissues with highest numbers and levels of lncRNA gene expression are the male reproductive organs. To learn why, we initiated a genome-wide analysis of Drosophila lncRNA spatial expression patterns in these tissues. The numbers of genes and levels of expression observed greatly exceed those previously reported, due largely to a preponderance of non-polyadenylated transcripts. In stark contrast to coding genes, the highest numbers of lncRNAs expressed are in post-meiotic spermatids. Correlations between expression levels, localization and previously performed genetic analyses indicate high levels of function and requirement. More focused analyses indicate that lncRNAs play major roles in evolution by controlling transposable element activities, Y chromosome gene expression and sperm construction. A new type of lncRNA-based particle found in seminal fluid may also contribute to reproductive outcomes.

Genome-Scale Investigation of the Regulation of azoR Expression in Escherichia coli Using Computational Analysis and Transposon Mutagenesis.

Bacterial azoreductases are enzymes that catalyze the reduction of ingested or industrial azo dyes. Although azoreductase genes have been well identified and characterized, the regulation of their expression has not been systematically investigated. To determine how different factors affect the expression of azoR, we extracted and analyzed transcriptional data from the Gene Expression Omnibus (GEO) resource, then confirmed computational predictions by quantitative reverse transcription polymerase chain reaction (qRT-PCR). Results showed that azoR expression was lower with higher glucose concentration, agitation speed, and incubation temperature, but higher at higher culture densities. Co-expression and clustering analysis indicated ten genes with similar expression patterns to azoR: melA, tpx, yhbW, yciK, fdnG, fpr, nfsA, nfsB, rutF, and chrR (yieF). In parallel, constructing a random transposon library in E. coli K-12 and screening 4320 of its colonies for altered methyl red (MR)-decolorizing activity identified another set of seven genes potentially involved in azoR regulation. Among these genes, arsC, relA, plsY, and trmM were confirmed as potential azoR regulators based on the phenotypic decolorization activity of their transposon mutants, and the expression of arsC and relA was confirmed, by qRT-PCR, to significantly increase in E. coli K-12 in response to different MR concentrations. Finally, the significant decrease in azoR transcription upon transposon insertion in arsC and relA (as compared to its expression in wild-type E. coli) suggests their probable involvement in azoR regulation. In conclusion, combining in silico analysis and random transposon mutagenesis suggested a set of potential regulators of azoR in E. coli.

Molecular Epidemiology of Carbapenem-Resistant Klebsiella aerogenes in Japan.

Information regarding Klebsiella aerogenes haboring carbapenemase in Japan is limited. A comprehensive nationwide survey was conducted from September 2014 to December 2022, and 67 non-duplicate strains of carbapenem-resistant K. aerogenes were isolated from 57 healthcare facilities in Japan. Through genetic testing and whole-genome sequencing, six strains were found to possess carbapenemases, including imipenemase (IMP)-1, IMP-6, New Delhi metallo-ß-lactamase (NDM)-1, and NDM-5. The strain harboring blaNDM-5 was the novel strain ST709, which belongs to the clonal complex of the predominant ST4 in China. The novel integron containing blaIMP-1 featured the oxacillinase-101 gene, which is a previously unreported structure, with an IncN4 plasmid type. However, integrons found in the strains possessing blaIMP-6, which were the most commonly identified, matched those reported domestically in Klebsiella pneumoniae, suggesting the prevalence of identical integrons. Transposons containing blaNDM are similar or identical to the transposon structure of K. aerogenes harboring blaNDM-5 previously reported in Japan, suggesting that the same type of transposon could have been transmitted to K. aerogenes in Japan. This investigation analyzed mobile genetic elements, such as integrons and transposons, to understand the spread of carbapenemases, highlighting the growing challenge of carbapenem-resistant Enterobacterales in Japan and underscoring the critical need for ongoing surveillance to control these pathogens.

The size diversity of the Pteridaceae family chloroplast genome is caused by overlong intergenic spacers.

BACKGROUND: While the size of chloroplast genomes (cpDNAs) is often influenced by the expansion and contraction of inverted repeat regions and the enrichment of repeats, it is the intergenic spacers (IGSs) that appear to play a pivotal role in determining the size of Pteridaceae cpDNAs. This provides an opportunity to delve into the evolution of chloroplast genomic structures of the Pteridaceae family. This study added five Pteridaceae species, comparing them with 36 published counterparts. RESULTS: Poor alignment in the non-coding regions of the Pteridaceae family was observed, and this was attributed to the widespread presence of overlong IGSs in Pteridaceae cpDNAs. These overlong IGSs were identified as a major factor influencing variations in cpDNA size. In comparison to non-expanded IGSs, overlong IGSs exhibited significantly higher GC content and were rich in repetitive sequences. Species divergence time estimations suggest that these overlong IGSs may have already existed during the early radiation of the Pteridaceae family. CONCLUSIONS: This study reveals new insights into the genetic variation, evolutionary history, and dynamic changes in the cpDNA structure of the Pteridaceae family, providing a fundamental resource for further exploring its evolutionary research.

Investigating the impact of insertion sequences and transposons in the genomes of the most significant phytopathogenic bacteria.

Genetic variability in phytopathogens is one of the main problems encountered for effective plant disease control. This fact may be related to the presence of transposable elements (TEs), but little is known about their role in host genomes. Here, we performed the most comprehensive analysis of insertion sequences (ISs) and transposons (Tns) in the genomes of the most important bacterial plant pathogens. A total of 35 692 ISs and 71 transposons were identified in 270 complete genomes. The level of pathogen-host specialization was found to be a significant determinant of the element distribution among the species. Some Tns were identified as carrying virulence factors, such as genes encoding effector proteins of the type III secretion system and resistance genes for the antimicrobial streptomycin. Evidence for IS-mediated ectopic recombination was identified in Xanthomonas genomes. Moreover, we found that IS elements tend to be inserted in regions near virulence and fitness genes, such ISs disrupting avirulence genes in X. oryzae genomes. In addition, transcriptome analysis under different stress conditions revealed differences in the expression of genes encoding transposases in the Ralstonia solanacearum, X. oryzae, and P. syringae species. Lastly, we also investigated the role of Tns in regulation via small noncoding regulatory RNAs and found these elements may target plant-cell transcriptional activators. Taken together, the results indicate that TEs may have a fundamental role in variability and virulence in plant pathogenic bacteria.

Remarkably High Repeat Content in the Genomes of Sparrows: The Importance of Genome Assembly Completeness for Transposable Element Discovery.

Transposable elements (TE) play critical roles in shaping genome evolution. Highly repetitive TE sequences are also a major source of assembly gaps making it difficult to fully understand the impact of these elements on host genomes. The increased capacity of long-read sequencing technologies to span highly repetitive regions promises to provide new insights into patterns of TE activity across diverse taxa. Here we report the generation of highly contiguous reference genomes using PacBio long-read and Omni-C technologies for three species of Passerellidae sparrow. We compared these assemblies to three chromosome-level sparrow assemblies and nine other sparrow assemblies generated using a variety of short- and long-read technologies. All long-read based assemblies were longer (range: 1.12 to 1.41 Gb) than short-read assemblies (0.91 to 1.08 Gb) and assembly length was strongly correlated with the amount of repeat content. Repeat content for Bell's sparrow (31.2% of genome) was the highest level ever reported within the order Passeriformes, which comprises over half of avian diversity. The highest levels of repeat content (79.2% to 93.7%) were found on the W chromosome relative to other regions of the genome. Finally, we show that proliferation of different TE classes varied even among species with similar levels of repeat content. These patterns support a dynamic model of TE expansion and contraction even in a clade where TEs were once thought to be fairly depauperate and static. Our work highlights how the resolution of difficult-to-assemble regions of the genome with new sequencing technologies promises to transform our understanding of avian genome evolution.

The evolution of transposable elements in Brachypodium distachyon is governed by purifying selection, while neutral and adaptive processes play a minor role.

Understanding how plants adapt to changing environments and the potential contribution of transposable elements (TEs) to this process is a key question in evolutionary genomics. While TEs have recently been put forward as active players in the context of adaptation, few studies have thoroughly investigated their precise role in plant evolution. Here, we used the wild Mediterranean grass Brachypodium distachyon as a model species to identify and quantify the forces acting on TEs during the adaptation of this species to various conditions, across its entire geographic range. Using sequencing data from more than 320 natural B. distachyon accessions and a suite of population genomics approaches, we reveal that putatively adaptive TE polymorphisms are rare in wild B. distachyon populations. After accounting for changes in past TE activity, we show that only a small proportion of TE polymorphisms evolved neutrally (<10%), while the vast majority of them are under moderate purifying selection regardless of their distance to genes. TE polymorphisms should not be ignored when conducting evolutionary studies, as they can be linked to adaptation. However, our study clearly shows that while they have a large potential to cause phenotypic variation in B. distachyon, they are not favored during evolution and adaptation over other types of mutations (such as point mutations) in this species.

Pangenome analysis reveals transposon-driven genome evolution in cotton.

BACKGROUND: Transposable elements (TEs) have a profound influence on the trajectory of plant evolution, driving genome expansion and catalyzing phenotypic diversification. The pangenome, a comprehensive genetic pool encompassing all variations within a species, serves as an invaluable tool, unaffected by the confounding factors of intraspecific diversity. This allows for a more nuanced exploration of plant TE evolution. RESULTS: Here, we constructed a pangenome for diploid A-genome cotton using 344 accessions from representative geographical regions, including 223 from China as the main component. We found 511 Mb of non-reference sequences (NRSs) and revealed the presence of 5479 previously undiscovered protein-coding genes. Our comprehensive approach enabled us to decipher the genetic underpinnings of the distinct geographic distributions of cotton. Notably, we identified 3301 presence-absence variations (PAVs) that are closely tied to gene expression patterns within the pangenome, among which 2342 novel expression quantitative trait loci (eQTLs) were found residing in NRSs. Our investigation also unveiled contrasting patterns of transposon proliferation between diploid and tetraploid cotton, with long terminal repeat (LTR) retrotransposons exhibiting a synchronized surge in polyploids. Furthermore, the invasion of LTR retrotransposons from the A subgenome to the D subgenome triggered a substantial expansion of the latter following polyploidization. In addition, we found that TE insertions were responsible for the loss of 36.2% of species-specific genes, as well as the generation of entirely new species-specific genes. CONCLUSIONS: Our pangenome analyses provide new insights into cotton genomics and subgenome dynamics after polyploidization and demonstrate the power of pangenome approaches for elucidating transposon impacts and genome evolution.

Identification and characterization of endogenous retroviruses upon SARS-CoV-2 infection.

Endogenous retroviruses (ERVs) derived from the long terminal repeat (LTR) family of transposons constitute a significant portion of the mammalian genome, with origins tracing back to ancient viral infections. Despite comprising approximately 8% of the human genome, the specific role of ERVs in the pathogenesis of COVID-19 remains unclear. In this study, we conducted a genome-wide identification of ERVs in human peripheral blood mononuclear cells (hPBMCs) and primary lung epithelial cells from monkeys and mice, both infected and uninfected with SARS-CoV-2. We identified 405, 283, and 206 significantly up-regulated transposable elements (TEs) in hPBMCs, monkeys, and mice, respectively. This included 254, 119, 68, and 28 ERVs found in hPBMCs from severe and mild COVID-19 patients, monkeys, and transgenic mice expressing the human ACE2 receptor (hACE2) and infected with SARS-CoV-2. Furthermore, analysis using the Genomic Regions Enrichment of Annotations Tool (GREAT) revealed certain parental genomic sequences of these up-regulated ERVs in COVID-19 patients may be involved in various biological processes, including histone modification and viral replication. Of particular interest, we identified 210 ERVs specifically up-regulated in the severe COVID-19 group. The genes associated with these differentially expressed ERVs were enriched in processes such as immune response activation and histone modification. HERV1_I-int: ERV1:LTR and LTR7Y: ERV1:LTR were highlighted as potential biomarkers for evaluating the severity of COVID-19. Additionally, validation of our findings using RT-qPCR in Bone Marrow-Derived Macrophages (BMDMs) from mice infected by HSV-1 and VSV provided further support to our results. This study offers insights into the expression patterns and potential roles of ERVs following viral infection, providing a valuable resource for future studies on ERVs and their interaction with SARS-CoV-2.

A concerted increase in readthrough and intron retention drives transposon expression during aging and senescence.

Aging and senescence are characterized by pervasive transcriptional dysfunction, including increased expression of transposons and introns. Our aim was to elucidate mechanisms behind this increased expression. Most transposons are found within genes and introns, with a large minority being close to genes. This raises the possibility that transcriptional readthrough and intron retention are responsible for age-related changes in transposon expression rather than expression of autonomous transposons. To test this, we compiled public RNA-seq datasets from aged human fibroblasts, replicative and drug-induced senescence in human cells, and RNA-seq from aging mice and senescent mouse cells. Indeed, our reanalysis revealed a correlation between transposons expression, intron retention, and transcriptional readthrough across samples and within samples. Both intron retention and readthrough increased with aging or cellular senescence and these transcriptional defects were more pronounced in human samples as compared to those of mice. In support of a causal connection between readthrough and transposon expression, analysis of models showing induced transcriptional readthrough confirmed that they also show elevated transposon expression. Taken together, our data suggest that elevated transposon reads during aging seen in various RNA-seq dataset are concomitant with multiple transcriptional defects. Intron retention and transcriptional readthrough are the most likely explanation for the expression of transposable elements that lack a functional promoter.

Cytogenomic Characterization of Transposable Elements and Satellite DNA in Passiflora L. Species.

The species Passiflora alata, P. cincinnata, and P. edulis have great economic value due to the use of their fruits for human consumption. In this study, we compared the repetitive genome fractions of these three species. The compositions of the repetitive DNA of these three species' genomes were analyzed using clustering and identification of the repetitive sequences with RepeatExplorer. It was found that repetitive DNA content represents 74.70%, 66.86%, and 62.24% of the genome of P. alata, P. edulis, and P. cincinnata, respectively. LTR Ty3/Gypsy retrotransposons represent the highest genome proportions in P. alata and P. edulis, while Ty1/Copia comprises the largest proportion of P. cincinnata genome. Chromosomal mapping by Fluorescent In Situ Hybridization (FISH) showed that LTR retrotransposons have a dispersed distribution along chromosomes. The subtelomeric region of chromosomes is where 145 bp satellite DNA is located, suggesting that these elements may play important roles in genome structure and organization in these species. In this work, we obtained the first global characterization of the composition of repetitive DNA in Passiflora, showing that an increase in genome size is related to an increase in repetitive DNA, which represents an important evolutionary route for these species.

Determinants of mer Promoter Activity from Pseudomonas aeruginosa.

Since the MerR family is known for its special regulatory mechanism, we aimed to explore which factors determine the expression activity of the mer promoter. The Tn501/Tn21 mer promoter contains an abnormally long spacer (19 bp) between the -35 and -10 elements, which is essential for the unique DNA distortion mechanism. To further understand the role of base sequences in the mer promoter spacer, this study systematically engineered a series of mutant derivatives and used luminescent and fluorescent reporter genes to investigate the expression activity of these derivatives. The results reveal that the expression activity of the mer promoter is synergistically modulated by the spacer length (17 bp is optimal) and the region upstream of -10 (especially -13G). The spacing is regulated by MerR transcription factors through symmetrical sequences, and -13G presumably functions through interaction with the RNA polymerase sigma-70 subunit.

Flexible TAM requirement of TnpB enables efficient single-nucleotide editing with expanded targeting scope.

TnpBs encoded by the IS200/IS605 family transposon are among the most abundant prokaryotic proteins from which type V CRISPR-Cas nucleases may have evolved. Since bacterial TnpBs can be programmed for RNA-guided dsDNA cleavage in the presence of a transposon-adjacent motif (TAM), these nucleases hold immense promise for genome editing. However, the activity and targeting specificity of TnpB in homology-directed gene editing remain unknown. Here we report that a thermophilic archaeal TnpB enables efficient gene editing in the natural host. Interestingly, the TnpB has different TAM requirements for eliciting cell death and for facilitating gene editing. By systematically characterizing TAM variants, we reveal that the TnpB recognizes a broad range of TAM sequences for gene editing including those that do not elicit apparent cell death. Importantly, TnpB shows a very high targeting specificity on targets flanked by a weak TAM. Taking advantage of this feature, we successfully leverage TnpB for efficient single-nucleotide editing with templated repair. The use of different weak TAM sequences not only facilitates more flexible gene editing with increased cell survival, but also greatly expands targeting scopes, and this strategy is probably applicable to diverse CRISPR-Cas systems.

Emergence of enhancers at late DNA replicating regions.

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

Earl Grey: A Fully Automated User-Friendly Transposable Element Annotation and Analysis Pipeline.

Transposable elements (TEs) are major components of eukaryotic genomes and are implicated in a range of evolutionary processes. Yet, TE annotation and characterization remain challenging, particularly for nonspecialists, since existing pipelines are typically complicated to install, run, and extract data from. Current methods of automated TE annotation are also subject to issues that reduce overall quality, particularly (i) fragmented and overlapping TE annotations, leading to erroneous estimates of TE count and coverage, and (ii) repeat models represented by short sections of total TE length, with poor capture of 5' and 3' ends. To address these issues, we present Earl Grey, a fully automated TE annotation pipeline designed for user-friendly curation and annotation of TEs in eukaryotic genome assemblies. Using nine simulated genomes and an annotation of Drosophila melanogaster, we show that Earl Grey outperforms current widely used TE annotation methodologies in ameliorating the issues mentioned above while scoring highly in benchmarking for TE annotation and classification and being robust across genomic contexts. Earl Grey provides a comprehensive and fully automated TE annotation toolkit that provides researchers with paper-ready summary figures and outputs in standard formats compatible with other bioinformatics tools. Earl Grey has a modular format, with great scope for the inclusion of additional modules focused on further quality control and tailored analyses in future releases.

A metagenomics pipeline reveals insertion sequence-driven evolution of the microbiota.

Insertion sequence (IS) elements are mobile genetic elements in bacterial genomes that support adaptation. We developed a database of IS elements coupled to a computational pipeline that identifies IS element insertions in the microbiota. We discovered that diverse IS elements insert into the genomes of intestinal bacteria regardless of human host lifestyle. These insertions target bacterial accessory genes that aid in their adaptation to unique environmental conditions. Using IS expansion in Bacteroides, we show that IS activity leads to the insertion of "hot spots" in accessory genes. We show that IS insertions are stable and can be transferred between humans. Extreme environmental perturbations force IS elements to fall out of the microbiota, and many fail to rebound following homeostasis. Our work shows that IS elements drive bacterial genome diversification within the microbiota and establishes a framework for understanding how strain-level variation within the microbiota impacts human health.

The curious case of the disappearing piRNAs.

Small non-coding RNAs are key regulators of gene expression across eukaryotes. Piwi-interacting small RNAs (piRNAs) are a specific type of small non-coding RNAs, conserved across animals, which are best known as regulators of genome stability through their ability to target transposable elements for silencing. Despite the near ubiquitous presence of piRNAs in animal lineages, there are some examples where the piRNA pathway has been lost completely, most dramatically in nematodes where loss has occurred in at least four independent lineages. In this perspective I will provide an evaluation of the presence of piRNAs across animals, explaining how it is known that piRNAs are missing from certain organisms. I will then consider possible explanations for why the piRNA pathway might have been lost and evaluate the evidence in favor of each possible mechanism. While it is still impossible to provide definitive answers, these theories will prompt further investigations into why such a highly conserved pathway can nevertheless become dispensable in certain lineages. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.